UCT-TP-262/04 
September 2004 



Electromagnetic proton form factors in large N c QCD* 

C. A. Dominguez, T. Thapedi 

Institute of Theoretical Physics and Astrophysics 
University of Cape Town, Rondebosch 7700, South Africa 

Abstract 

The electromagnetic form factors of the proton are obtained using a particular realization of QCD in 
the large N c limit (QCD^), which sums up the infinite number of zero- width resonances to yield an 
Euler's Beta function (Dual-QCD^). The form factors F\(q 2 ) and F2(q 2 ), as well as Gm(? 2 ) agree 
very well with reanalyzed space-like data in the whole range of momentum transfer. In addition, the 
predicted ratio [i v Ge/Gm is in good agreement with recent polarization transfer measurements at 
Jefferson Lab. 
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Quantum Chromodynamics (QCD) in the limit of large number of colours (QCD^) [1] is known to predict 
a hadronic spectrum consisting of an infinite number of zero- width resonances [2]. However, since real QCD 
has never been solved exactly and analytically, the hadronic parameters (masses, couplings, etc.) remain 
unpredicted. A few models of this spectrum have been proposed for heavy quark Green's functions [3]- [4], 
as well as for light quark systems [5]. The infinite number of zero-width resonances of QCD^ is reminis- 
cent of Veneziano's dual-resonance model [6], the precursor of string theory. In fact, inspiration from this 
model has led to a proposal called Dual-QCD^ [7], a specific realization of QCD^ where the masses and 
couplings in a Green's function are chosen to yield an Euler's Beta function of the Vencziano type. For 
three-point functions, the form factors exhibit asymptotic Regge-behaviour, i.e. power-behaviour, in the 
space-like region controlled by a single free parameter. Dual-QCD^ has been applied quite successfully 
to the electromagnetic form factor of the pion in the space-like region [7]. Indeed, results are in excellent 
agreement with experiment, far better than e.g. naive Vector Meson (rho-) Dominance or purely pertur- 
bative QCD [8]. This is the case not only for the pion form factor itself, but also for the mean-square 
radius, and the observed deviation from universality (the ratio g p7T7r / ' f p ) . In addition, unitarization of 
Dual-QCD^ leads to a prediction of the vector two-point spectral function, in the time-like region, in 
reasonable agreement with data (a more refined model in the time-like region has been proposed recently 
in [9]). Encouraged by this success, we discuss in this note an analysis of the electromagnetic proton form 
factors in the framework of Dual-QCD^. 



The Dirac and Pauli form factors of the proton, F 1 (q 2 ) and F 2 (q 2 ), respectively, are defined as 

< N(p 2 )\J* M (0)\N( Pl ) >= u N {p 2 ) [F^h, + ^-F 2 {q 2 )a^] u N ( Pl ) , (1) 

where q 2 = (p 2 —pi) 2 , while k = fi p — l, and M v are the proton's magnetic moment and mass, respectively, 
and Fi i2 (q 2 ) are normalized as Fi(0) = F 2 (0) = 1. On the other hand, the Sachs' form factors Gs(</ 2 ) and 
G M {q 2 ) are given by 

G E {q 2 ) = Fi(q 2 ) — k t F 2 (q 2 ) 

G M (q 2 ) = Fi(q 2 ) + k F 2 (q 2 ) , [Z > 

where r = -q 2 /4M 2 = Q 2 /AM 2 , and the normalization is then G E (0) = 1, and G M (0) = [i v - 
In the very early applications of the dual-resonance model to three point functions involving more than 
one form factor [10], it was not quite clear to which form factor, or linear combination of form factors, 
should the model apply. There is no ambiguity in Dual-QCD^, as this is a realization of a quantum field 
theory. The form factors should then be those appearing in the primary hadronic spectral function, dual 
to the QCD field theory spectral function. In other words, the form factors with the correct pole structure 
satisfying dispersion relations in the complex energy plane. In the case of the nucleon, these are the Dirac 
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and Pauli form factors which in QCD^ become 



n=0 

where s = q 2 , and the masses of the vector- meson zero- width resonances, M n , as well as their couplings 
C\ n and C 2 „, are not predicted. In Dual-QCD^ these are chosen so that the form factors become Beta 
functions (ratios of gamma- functions) , i.e. 

_ r(/3 1|2 - 1/2) (-1)" i u , 

(1,2)n r(n + l) r(/3 1>2 -l-n) ' [ > 

where (3\ t 2 are free parameters controlling the asymptotic behaviour in the space-like region (s < 0), and 
a' = l/2M p 2 is the universal string tension in the rho-meson Regge trajectory 

a p (s) = 1 + a'(s - M 2 ) . (5) 

The mass spectrum is chosen as [11] 

M 2 = M 2 (l + 2n). (6) 
Using Eqs.(4) and (6) in Eq.(3) one obtains 

rC9i, 2 -l/2) ^ i i 



F iAs) - -j= X. 



where B(x,y) is Euler's Beta function. In the time- like region (s > 0) the poles of the Beta function 
correspond to an infinite set of zero- width resonances with equally spaced squared masses given by Eq.(5). 
In fact, from Eq.(7) it follows 

Im FlM ■ m^i g Mr i S(Mi _ . (8) 

v n— ' 

Asymptotically, the form factors in the space-like region exhibit Regge-behaviour, viz. 

lim F 12 (s) = (-a' s )( 1 -/ ? i^) , (9) 

The free parameters /?i ;2 can be fixed from fits to the data in the space-like region. Notice that the values 
(3i t 2 = 2 reduce the form factors to single rho-meson dominance (naive Vector Meson Dominance). The 
mass formula Eq.(6) predicts e.g. for the first three radial excitations: M p > ~ 1340 MeV, M p » ~ 1720 
MeV, and M p ,„ ~ 2034 MeV in reasonable agreement with experiment [12] : M p > = 1465 ± 25 MeV, 
Mp» = 1700 ± 20 MeV, and M p », = 2149 ± 17 MeV. Alternative (non-linear) mass formulas might be 
required if one were to match the asymptotic Regge behaviour to the Operator Product Expansion of 
current correlators at short distances [13]. However, the differences in the values of the first few masses 
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are at the level of a few percent. Hence, the form factors would hardly be affected, since the contribution 
from high mass states is factorially suppressed. 



Historically, the Sachs form factors were first determined from measurements of elastic electron-proton 
scattering cross sections (Roscnbluth technique) [14]. Direct extractions of G_e(q 2 ) and Gm{<1 2 ) up to 
—q 2 = Q 2 ~ 7 GeV 2 indicated the empirical approximate scaling relation: h p Ge{(1 2 )/Gm{(1 2 ) — 1- At 
higher values of Q 2 , the contribution of Gb((? 2 ) to the cross section is kinematically suppressed. On the 
other hand, recent electron-proton polarization transfer measurements at Jefferson Lab (JLab) [15] up to 
Q 2 ~ 6 GeV 2 show a considerable deviation from this scaling behaviour, except possibly at very small 
Q 2 [16]. After some debate about the source of the discrepancy between cross-section (Roscnbluth) and 
polarization transfer extractions of the form factors [17], it appears that the culprit is the two-photon 
exchange correction [18]. We assume this to be the case, and adopt the experimental data on the Sachs 
form factors as corrected in [19]. These corrections are made in order to bring the Rosenbluth data into 
agreement with the polarization transfer data on the ratio ^ p Gb((7 2 )/G'm((? 2 )- We then use Eq.(2) to 
obtain data points for Fi(q 2 ) and F 2 (q 2 ). After fitting this data base with Eq.(7) we find 0i = 3.03, and 
02 = 4.20. Figures 1, and 2 show the results of the fits for i*\(g 2 ), and F 2 (q 2 ), corresponding to these values 
of 0i,2, together with the corrected data points of [19]. In Fig. 3 we show Gm(<? 2 ), as obtained from Eq.(2) 
using the fitted Fi y2 (q 2 ) 7 together with the same data base. As can be appreciated, the agreement between 
Eq.(7) and the data is very good. Having fitted a data base corrected to account for the polarization 
transfer data on the ratio fj, p GE(q 2 )/GM(q 2 ), we would expect our theoretical form factors to lead to a 
ratio in agreement with experiment. While this is the case, it turns out that this ratio is very sensitive 
to the pair of values 0i — 2 , with a strong correlation between them. In Fig. 4 we show the theoretical 
prediction of the ratio /x p G_e(9 2 )/Gm(<7 2 ) corresponding to 0i = 3.0 and 02 — 4.2, together with the JLab 
data [15]. Small variations of these parameters lead to correlated pairs resulting in equally good fits, e.g. 
the pair 0i = 2.95 and 2 = 4.13 leads to an almost identical theoretical prediction. Exploring the 0i-0 2 
parameter space, and performing a combined fit to Fi(q 2 ), F 2 (q 2 ), and the ratio n P GE(q 2 )/GM(q 2 ) gives 



where ip(x) is the digamma function. Using the results from Eq.(10) in Eq. (11) gives < r\ > x / 2 — 0.72 fm, 
and < r\ > 1 / 2 = 0.78 fm. From Eq.(2), these radii lead to the Sachs radii < r| > 1 / 2 = 0.81 fm, and 
< r 2 M > 1 / 2 = 0.76 fm. These values are in reasonable agreement with various other results in the literature 



(3i = 2.95-3.03 
02 = 4.13 - 4.20 . 



(10) 



The mean-squared electromagnetic radii that follow from Eq.(7) are given by 



< r\ 2 >= 6 a' [V(/?i,2 - 1/2) - ^(1/2)] , 



(11) 
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[20], especially taking into account that the free parameters j3\^ have been fixed from the large-Q 2 data, 
while the form factors decrease by 3-4 orders of magnitude in the range < Q 2 ?S 30 GeV 2 . It should also 
be kept in mind that the strong deviation from unity of the ratio [i v G e(q 2 ) / 'G m(<? 2 ) might affect existing 
extractions of the Sachs radii from data. 

In summary, the nuclcon form factors Fi(q 2 ), and F 2 (q 2 ) obtained in the framework of Dual-QCD^ 
reproduce very nicely the experimental data in the space-like region, as corrected in [19], in the wide 
range < Q 2 ^ 30 GeV 2 . The Sachs magnetic form factor, Gm(? 2 ), as well as the non-trivial ratio 
HpGE{q 2 )/GM{q 2 ) can also be accounted for in this framework. These results provide strong support for 
Dual-QCD^ as a viable realization of QCD in the large N c limit. 
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Figure Captions 

Figure 1. Dual-QCDoo form factor Fi(Q 2 ), Eq.(7), for the fitted parameter f}\ — 3.03, together with the 
experimental data as corrected in [19]. 

Figure 2. Dual-QCDoo form factor F 2 (Q 2 ), Eq.(7), for the fitted parameter /3 2 = 4.20, together with the 
experimental data as corrected in [19]. 

Figure 3. Dual-QCDoo form factor Gm(Q 2 ), Eq.(2), for the fitted parameters (3\ — 3.03, and (32 — 4.20, 
together with the experimental data as corrected in [19]. 

Figure 4. Dual-QCDoo ratio h p G e (Q 2 )/Gm(Q 2 ) for the fitted parameters /?i = 3.00, and (3 2 = 4.20, 
together with the experimental data [15]. 
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Figure 4: 
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